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The influence of fallback discs on the spectral and timing 
properties of neutron stars 
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ABSTRACT 

Fallback discs around neutron stars (NSs) are believed to be an expected outcome of 
supernova explosions. Here we investigate the consequences of such a common out- 
come for the timing and spectral properties of the associated NS population, using 
Monte Carlo population synthesis models. We find that the long-term torque exterted 
by the fallback disc can substantially influence the late-time period distribution, but 
with quantitative differences which depend on whether the initial spin distribution is 
dominated by slow or fast pulsars. For the latter, a single-peaked initial spin distribu- 
tion becomes bimodal at later times. Timing ages tend to understimate the real age of 
older pulsars, and overestimate the age of younger ones. Braking indices cluster in the 



range 1.5 



3 for slow-born pulsars, and —0.5 



5 for fast-born pulsars, with 



the younger objects found predominantly below n ^ 3. Large values of n, while not 
common, are possible, and associated with torque transitions in the NS+disc system. 
The 0.1-10 keV thermal luminosity of the NS+disc system is found to be generally 
dominated by the disc emission at early times, t < 10 3 yr, but this declines faster than 
the thermal surface emission of the neutron star. Depending on the initial parameters, 
there can be occasional periods in which some NSs switch from the propeller to the 
accretion phase, increasing their luminosity up to the Eddington limit for ~ 10 3 -10 4 
years. 

Key words: accretion, accretion discs - pulsars: general. 



1 INTRODUCTION 

Neutron stars (NSs) are born in supernova (SN) explosions. 
While the stellar core collapses, the envelope is ejected out- 
wards. The ejected material carries angular momentum, due 
■ to the rotation of the progenitor star. If a fraction of the 
ejecta does not have enough energy to escape, it w ill form 
a fallback accretion d isc around the collapsed core (|Michell 

I198SI ; IChevalierlll989h . Although a systematic and quanti- 
tative study of the fate of the ejected material following 
a SN explosion is still lacking, fallback discs around NSs 
have been invoked to explain a wide range of NS properties, 
such as enh anced X-ray luminosit i es in the Anomalous X- 
ray Pulsar ijCh attcrjcc ct al. 2000; Ichatteriee fc Hernquistl 

l200d ; I AlparlboOll ; lAlpar et al.ll200llh pulsar braking i ndices 
that suggest non-dipolar torques ( Menou et al .112001 al l, pul- 
sar jets (Blackman & Perna 2004; Blaschke, Grigorian & 
Voskresensky 2004), discrepancies between the real (kine- 



matic) pulsar ages an d their characteristic ages r = P/2P 
jMarsden et al.ll2001al lbl). 

From an observational point of view, the search for 
fallback discs has been quite active in the last decade 
or so. Following up on the suggestion that fallback discs 
around NSs are brighter in the infrared (IR) and at longe r 
wavelengths (| Perna fc Hernquistl |2000| ; iPerna et al.l 120001 ). 
a number of isolated NSs hav e been imaged in tho se 



Kaplan et alJl200ll; J l srael et alj|2003l ; llvlignani et alj 12007 



Wan g et al 



Wang et al 



K US 

2007al lbJ; IPosselt et all l201Ch . In particular, 
( 20061 1 reported the first direct detection of a 
fallback disc around an isolated NS, based on an IR excess 
which could be well modelled with emission from a passive, 
dusty disc. Eviden ce for IR excess has also been reported in 
other NS systems jlsrael et al.ll2003l ; [Kaplan et alj|2009ft . 

If fallback discs around NSs are indeed a ubiquitous 
phenomenon, they can substantially interfere with the spin 
evolution of the NS. This will affect estimates of the NS 
spin birth distribution based on the measured values of the 
period P and its derivative P, and the simple assumption 



bands dCoe fc PightlinS Il998); iHulleman et all l2000a 
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that the pulsar slow down is due to dipole radiation only. 
The distribution of NS spin births, which contains informa- 
tion on the physics of the SN explosion, has indeed been 
the subject of of e xt ensive in vestigation for several d ecades 



[1973; 



Phinnev t 
19851; 



Blandfordl Il98ll: 



llGunn fc Ostrikcr 

Srinivasan et al.lll984l;lL~vne et al.ll 19851: iNaravan fc Ostriker 
1990; Arzoumanian et al.l 120021 ; iFaucher-Giguere fc Kaspil 
2006 ; [Perna et al.ll200 



In this paper we perform a statistical study of the spin 
evolution of NSs under the combined evolution resulting 
from dipole radiation losses and the interaction with a fall- 
back disc. We aim at a) assessing to what extent the inferred 
spin birth properties of pulsars can be influenced by the pres- 
ence of fallback discs, b) identifying observational properties 
that can help confirm/rule out a ubiquitous presence of fall- 
back discs and c) exploring whether fallback discs can help 
explain some currently puzzling pulsar observations. 

Our paper is organized as follows: the coupled evolu- 
tion of the disc and pulsar system is described in Section 
2.1. Evolutionary tracks for several specific combinations of 
initial conditions are presented and discussed in Section 2.2 
(timing properties) and Section 2.3 (luminosity). We present 
the statistical timing and luminosity properties of the pulsar 
population, derived via Monte Carlo simulations, in Section 
3.1, while in Section 3.2 we further discuss the luminosity 
properties of the NS+disc system, and its differences from 
other astrophysical sources. We summarize and conclude in 
Section 4, highlighting puzzling pulsar observations which 
can be accounted for by the presence of fallback discs. 



torque is applied at Ri c . In summary, the inner disc radius 
is located at Ri n — min [Rm, -Ric]- 

Inspection of Eq.(JT| shows that, at any given time, if 
Q < Q,K(Rin), the disc exterts a positive torque, which, if 
dominant over the dipole component, would make the pul- 
sar spin up and accrete. On the other hand, if Q > QuiRin), 
mate rial from the disc is stopp ed by the centrifugal bar- 
rier (|lllarionov fc Sunvaevl [l~975l l . and the disc contributes 
to furher spin down the pulsar. 

Formation of the fallback disc ensue as material from 
the stellar envelope is thrown out at large distances during 
the supernova explosion; the fraction of this material which 
remains bound, falls back, and circularizes onto Keplerian 
orbits corresponding to their original specific angular mo- 
mentum (Chevalier 1989; MacFadyen et al. 2001; see also 
Perna & MacFadyen 2010 for an exploration of different an- 
gular momenta distributions in collapsing stars). 

Once the initial disc has been formed, its subsequent 
evolution can be mod elled following the res ults of the nu- 
merical simulations of ICannizzo et all (|l990t ). They showed 
that the accretion rate in the debris disc follows the evolu- 
tion 



M(t) = 



M , 

M (t/t ) 



for < t < t , 
for t > to, 



(2) 



where Mo is the accretion rate during an init ial plateau 
phase lasting a timescale tn. lMenou et al l (l2001bl ) estimated 
the timescale of the transient accretion phase to be on the 
order of 



to 



R 1 ^ 



6.6 x lO-^Tce)- 1 ^ 2 ^) 



yr, 



(3) 



2 A MODEL OF FALLBACK ACCRETION 

2.1 Time evolution of the accretion flow 

We model the combined dipole/disc torque acting on the 
N S (rotating at angul ar velocit y Q.) following the ap proach 
of iMenou et all l| 19991 ) (see also lMenou et alj|2001al ): 



IQ. = -/3fi 3 + 2MRf n n K (Riv 



1 - 



il K {Rin) 



(1) 



where Ri n is the inner radius of the disc, / = |Mns-Rns i s 
the moment of inertia of the NS, f3 = B 2 Rf iS /6c 3 (assum- 
ing for simplicity an orthogonal rotator), M is the accretion 
rate and fix is the Keplerian angular velocity, both evalu- 
ated at Ri n . The first term in the r.h.s. of Eq.|TJ) represents 
the torque due to dipole radiation losses, whereas the second 
term describes the torque exerted by the disc. At high accre- 
tion rates, the ram pressure of the accreting material causes 
the disc to penetrate inside the magnetosphere, up to an in- 
ner radius at which the magnetic pressure balances the ram 
pressure. This special location is called the magnetospheric 
radius, R m = 2.55 x 10 8 M~ 6 2/7 M~^ ' B% 7 cm. 

If, at any point during the coupled NS+disc evolution, 
the magnetospheric radius, i? m , formally exceeds the light 
cylinder radius, R\ c = c/Q, the field lines on the disc are 
no longer able to exert a significant pressure on the disc. 
On the other hand, the ram pressure of the material from 
the accreting disc continues to act. The disc is not able to 
penetrate inside R\ c , because the magnetic pressure from 
the NS would then push it outwards. Therefore, when J? m 
becomes formally larger than Ri c , we assume that the disc 



where c s is the sound speed, a' the viscosity parameter, 
Rd,s(to) the initial disc radius in units of 10 s cm, and T Cy0 
a typical disc temperature during the early phases, in units 
of 10 6 K. The index a of the powerlaw decay, was found 
to be fs 19/16. A similar behaviour, i.e. a plateau followed 
by a powerlaw decline, w as also found in the simulations by 
iMacFadven et al.l l|200ll) . although with a ~ 5/3. Here, we 
adopt a = 19/16 for our model. 

For an estimate of the expected range of values of the 
initial accretion rate, we rely on results of numerical sim- 
ulations of massive star collapse and subsequent fallback 
onto the compact object. Mineshige et al. (1997) found that 
the fallback material accretes towards the central object 
with a supercritical accretion rate, exceeding the Edding- 
ton limit by at least a factor of 10 6 . The simulations of 
IMacFadven et al.l l|200ll) . for the case of a 25 Mq star, yield 
Mo ~ 10 27 — 10 29 g s _1 . More recently, fallback in supernovae 
was studied by Zhang, Woosley & Heger (2008) for a wide 
range of masses, and different metallicities and explosion en- 
ergies. The full hydrodynamical study in one representative 
case yielded initial accretion rates of ~ 10 29 erg s -1 . An 
analysis of SN 1987A (Chevalier 1989) showed that the ac- 
cretion rate at the time of the reverse shock (7 x 10 3 s) is 
2.2 x 10 28 g s _1 . In general, when considering a distribution 
of NSs, a distribution of initial]] accretion rates is expected, 
as this depends on the angular momentum of the progenitor 



1 Note that here we speak of 'initial' disc accretion rate as if the 
disc had been formed instantaneously, since the timescale of disc 
formation, on the order of the dynamical timescale of collapse 
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star, its mass, and the energetics of the explosion. In our 
work, to be broadly consistent with the result of previous 
numerical work, as well as of previous semianalytical investi- 
gations which used fallback discs, we select the initial accre- 
tion rate Mo to be distributed log-uniformly from 10 25 g s _1 
to 10 28 g s _1 . However, with these values, the mass accre- 
tion inflow at early times would be highly supercritical; that 
is, it would produce a highly super-Eddington luminosity if 
all the inflowing matter was efficiently accreted through the 
disc and onto the NS. The classical Eddington limit for a 
NS, L Ed d = 1.3 x 10 38 (M NS /M ) 



2 x 10 38 erg s *, re- 
, for a radiatively 



-i 



10 1S gs 



quires only accretion rates M 
efficient flow. 

Faced with a highly supercritical accretion flow at the 
outer boundary, an ac creting system has only two options 
jPout anen et al.l 120071 ) : it can accrete most of it in a ra- 
diatively inefficient way (e.g., the "Polish donut" or "slim 
disc" classes of accretion flows: lAb ramowicz et all 1 19881 : 
IWatarai et all [2000 ) ; or it can spend most of the accretion 
power to blow awa y the excess infalling mass , in a radia- 
tively driven wind (Sh akura fc Sunvaevl 11973 ). For an ac- 
creting black hole, observational data and hydrodynamical 
simulations suggest th at the physical situation is a combi- 
nation of both effects ijOhsuea fc Mincshigel [201ll ) . For an 
accreting NS, the presence of a hard surface prevents advec- 
tive solutions: the excess gravitational energy released dur- 
ing accretion cannot be carried through a horizon and will 
always be dissipated, eventually. Therefore, in our model we 
assume that super-critical inf lows on a young NS will le ad to 
massive outflow s . Follo wing IShakura fe Sunvaevl l| 19731 ) and 
iPoutanen et al. (2007), we assume that only M ~ M^dd = 
iEdd/O.lc 2 reaches the inner edge of the accretion disc (from 
where it can then plunge onto the NS), regardless of the 
mass supply at the outer disc; and that the accretion rate 
through each radius decreases linearly, Al ~ R, because of 
radiation-driven wind losses. Thus, we parameterize the ra- 
dial dependence of the mass flow in the disc at a fixed time 
as 



M(R) = Al(R in ) + [M(Ro 



M(i? in 



j R — Rin 

Rout — R ! ) ) 



(4) 



where M(R out ,t) = M(R out ,t ), for < t < to, and 
Al (Rout, to) (t/to)~ a , for t > to, is the rate of fallback 
mass inflow at the outer boundary, and M(R ln ,t) = 

min j^A/Edd, Al(R ou t, f)j • In the rest of this paper, for clar- 
ity, we shall use the term "accretion rate" when referring 
to M(Ri n ) (what is actually accreted through the whole 
disc), and "inflow rate" for A1(R) at larger radii, most 
of which is lost in the outflow at early epochs. (See also 
iKuncic fc BicknelJl2007l for analytical modelling of accretion 
discs with a power-law scaling of the inflow rate M(R) ~ 
R b .) Recalling the definition of magnetospheric radius i? m , 
it is possible that at early times, when the accretion rate is 
near Eddington, i? m < _Rns : iu that case the disc extends all 
the way to the NS surface, and we assume the disc torque 

to be M E dd-RNS^K(-RNs)- 

At intermediate epochs, when M(R out ,t) w 



of the star envelope (tens of seconds), is much shorter than the 
timescale (thousands of years) of interest here for the fallback disc 
to affect the pulsar timing evolution. 
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Figure 1. Range of the P — B parameter space in which the 
pulsar is either spinning up or down, owing to the combined effect 
of dipole radiation losses and torque by a disc. For each value of 
the accretion rate at the magnetospheric radius, the region above 
its corresponding curve represents the parameter space where the 
pulsar is spun up, while the region below represents spin down 
conditions. 



M(Riu,t) < A^Edd, the wind subsides, and the inflow 
rate becomes constant thoughout the disc and equal to the 
accretion rate at the inner edge of the disc. Finally, at late 
times, the accretion rate declines further, until it reaches 
a value comparable to what is expected from accretion 
from the interstellar medium. From that epoch onwards, we 

equal to the 



take Al to be constant^ Al 



2 x 10 n _s_s _1 , 



Bondi accretion rate l|Bondilll952l : iFrank et al.ll2002T l for a 
medium with characteristic temperature kT ~ 10 3 keV and 
number density n ~ 1 cm -3 . Although the precise value of 
the Bondi accretion rate in a low de nsity medium depends 
on the details of the accretion mode ( Blaes fc Madaulll993l : 
IPopov et all l2000l : iPerna et ail 120031 ). our results, which 
focus on moderately young pulsars, are not very sensitive 
to the choice of the late-time asymptotic value. 



2.2 Time evolution of the pulsar timing properties 

As discussed above, at any time during the pulsar evolution, 
the disc can either contribute to spinning down the pulsar 
even more (with respect to the dipole term alone), or it can 
offset the dipole term and spin the pulsar up. Fig. 1 shows, 
for several values of the accretion rate, the P — B parameter 
space in which the pulsar spins up (down). 

At late times, the torque may disappear if the disc be- 
comes neutra l , viscou s dissipation drops and accretion stops. 
iMenou et al.l l|2001bh estimated a neutralization timescale 



tn — 



-Rd,lo(*o) 



M(t )/10 16 gs- 



to, 



(5) 



where Rd,io is the disc radius in units of 10 cm. The pos- 
sible onset of a transition to a phase of 'dead' disc is rather 
uncertain, as it depends on the nature of viscosity; opacity 
(and hence disc composition) also plays an important role, 
as well as the disc size (Menou et al. 2001b). Furthermore, 
it has been pointed out that irradiation by the pulsar is 
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likely to keep the disc ionized and continue the mass inflow 
(Alpar et al. 2001). Here, to avoid introducing highly un- 
certain parameterizations, we assume that the disc is active 
throughout the active lifetime of the pulsar. If the disc were 
to become dead at some point, from that time on the pulsar 
would simply continue its evolution as an isolated pulsar, 
spinning down by dipole radiation losses alone (neglecting 
gravitational energy losses and other secondary effects). 

The birth distribution of surface magnetic fields and 
spin periods has been the subject of numerous studies. Usu- 
ally, a log-Gaussian distribution is assumed for both. For the 
B field distribution, there is no w a general consensu s abou t 
the initial parameters; following I Arzoumanian et all < |2002h . 
we take 

(logSo(G)) = 12.35, a logSo = 0.4. (6) 

Whether accretion leads to field decay, and the timescale 
over which the magnetic field re-emerges, are still unsettled 
issues, especially with regards to quantitative details. Sev- 
eral papers have been devoted to this topic (e.g. Cumming et 
al. 2001; Geppert et al. 1999), and the coupled evolution be- 
tween a NS with a fallback disc and a field gradually buried 
by the disc accretion has been studied by Liu & Li (2009), 
using an emphirical expression for B-field decay presented 
by Shibanazi et al. (1989), B = B /(l + AM/m B ), where 
A M is the accreted mass, and tub is a mass normalization 
constant with typical values between ~ 10~ 5 and 10 _4 Mq. 
The resulting evolutionary tracks for the periods of the NS 
studied were found to be very model-dependent, with the 
magnitude of the effect dependent on both the initial B 
field and the disc properties. In the current work, to avoid 
results dependent on a rather uncertain parametrization, we 
consider the case of negligible field burying. We stress that 
our investigation is aimed at capturing features of pulsar 
properties which could be attributed to NS-disc interaction, 
but we do not make any attempt to quantitatively model 
data. 

For the birth period distribution, on the other hand, 
there is a widely different range of values inferred by different 
investigations: from a dist ribution peaking around a few ms 
|Arzoumanian et al.ll2002T ) 

(logPo(s)) = -2.3, <7 logPo > 0.2, (7) 

to one peaking at several hundred ms (Faucher-Giguere & 
Kaspi 2006) 

(logPo(s)) = -0.5, <j logFo = 0.2, (8) 

A lower limit of several tens of ms was derived from the 
scarcity of bright X-ray pulsa rs associated with historical 
supernovae |Perna et al.l I2008T ). It should be pointed out, 
however, that all those investigations assumed that pulsar 
evolution occurs with no contribution from fallback torques. 
Therefore, we use them here as an initial reference point to 
study the effect that hypothetical fallback discs would have 
on the pulsar evolution, at a statistical level. We take an 
initial spin distribution that encompasses the broad range 
of values inferred by previous investigations. In fact, one of 
the goals of our study is to uncover differences in the effect 
that a torque by a disc has on NSs born with fast and slow 
spins. However, we do need to note that it is possible that 
discs might be much more common around slow-born NSs 
than fast-born NSs, since the rotational energy of the latter 



can rapidly disrupt the disc (Eksi, Hernquist & Narayan 
2005). 

The presence of a fallback disc torque, in addition to 
the dipole contribution, can be tested with measurements of 
the braking index, an observable quantity defined as 

n = - — , (9) 

and of the timing age, 



where P = 2n/Q. If the primary source of spin-down is due 
to dipole radiation losses, then n = 3, and the timing age of 
the pulsar is equal to the real age, if Po <C P. 

Figure [2] shows the time evolution of the braking in- 
dex, timing age, and the absolute magnitude of the dipole 
and disc torques, for four representative pulsars, suitably 
chosen to represent a wide range of initial birth param- 
eters. For a fixed initial B field strength Bo (defined in 
Eq , we consider the time-evolution for a typical slow- 
born (Pq = 300 ms) and fast-born (Pq — 5 ms) pul- 
sar, with an initial fallback rate either at the maximum 
(Mo, max = 10 28 g s _1 ) or at the minimum value (Mo,min — 
10 25 g s" 1 ) of the assumed distribution. In order to com- 
pute the evolutionary tracks shown in Figure [2] by means of 
Eq.® and (|10p . we use Eq.(JT| to evaluate Cl, from which 
Q and Q are derived through numerical differentiation and 
numerical integration, respectively. For the latter, the ini- 
tial condition is fl(t = 0) = 2n/Po, with Po = Pq or 
Po = Po depending on the case. Note that in Eq.{T}, M 
represents the accretion rate at the inner disc boundary, 
i.e. M(Pi n ) = min ^M E dd, M (R ou t, t)j , where, recalling from 

Section 2.1, M(R out ,t) = M(R out ,t ), for < t < to and 
M (Pout, to) (t/to)~ a for t > to. For each combination of sys- 
tem parameters Po, Po, M (Rout, to), the timing evolution of 
the NS driven by the coupled disc+dipole torques is there- 
fore uniquely determined. 

The result of our calculations is that the disc torque is 
always negative, that is f2 > £7k (corresponding to a spin 
down/propeller regime), except for a few brief phases when 
the accretion rate is very high and £1 < Ok (correspond- 
ing to a spin up/accretor regime, see also Section 2.3). The 
absolute magnitudes of the disc and dipole torques (dotted 
and dashed lines, respectively) are plotted in Figure [2] the 
short epochs in which the sign of the disc torque is positive 
are plotted with red dots (Figure [2] left-hand panels). The 
presence of significant (and variable) disc torques causes the 
braking index to vary. More generally, n goes through dis- 
continuities caused by transitions between disc-driven spin- 
up and spin-down phases, between phases of constant accre- 
tion rate and powerlaw decay, and when the location of the 
inner radius of the disc switches between P m and P; c . In be- 
tween those rapid transitions, the complex behaviour of the 
braking index is mostly determined by the non-linear cou- 
pling between NS spin, inner radius of the disc and accretion 
rate. In particular, for fast pulsars, the dipole torque largely 
dominates over the disc torque at early times, while the lat- 
ter dominates at later times. Instead, for slow pulsars, the 
disc torque is larger than the dipole component at almost all 
times. Both slow and fast pulsars have a roughly constant 
braking index during the first few 10 3 yrs, corresponding to 



© 2012 RAS, MNRAS 000, [TV?? 



Fallback discs in neutron stars 5 



1 0' 




Pq • "»0,MAX 




cn 9 

s 1CT 2 






o 

— icr 4 

- 1CT 6 


"""----V- 

- Disk 

Dipole 




6 


: I — I — I — | I — t- 

f / 


j4-™| | i mill 

E 


A 

ZL 


\ J 






/ 




10 io 


i i — i — i — i — t 


L™i — I — ] 






1 s 






^3 10 6 

10 4 






1 o 2 







1 ! 




10°r 






10~ 2 r 


in 






1(r 4 r 


O 






10- 6 r 




10 -io r 



0° 10 1 



1 3 1 4 1 5 1 6 



o 7 



t (yr) 




10 6 10 7 






10 4 r 




1 2 r 


CO 


10°r 


cn 


10~ 2 r 


O 


1CT 4 r 




1CT 6 r 




10~ 8 r 




10 -io r 




Figure 2. Time evolution of the dipole and disc components of the spin-down rate f2 , of the braking index n, and of the timing age r, for 

slow-born (Pq = 300 ms; top panels) and fast-born (Pq = 5 ms; bottom panels) pulsars, starting with either a high (Mo, MAX = 10 28 g s _1 ; 
left-hand panels), or a low (Mq MIN = 10 25 g s — 1 ; right-hand panels) fallback inflow rate at the outer boundary of the disc. f2 is always 
negative (spin-down), except for occasional brief disc spin- up phases, marked by red dots (left-hand panels). For the breaking index, we 
note a characteristic n m 1 transitional phase that happens only in fast-born pulsars (bottom panels). 
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the epoch when M(R out ) w constant and A/(i?i n ) w A^Edd- 
However, the index is exactly n = 3 for (dipole-dominated) 
fast pulsars, while it is n ~ Q/[fi — fix (-Rin)] for (disc- 
dominated) slow pulsars (from Eq.(l)). The constancy of n 
in the regime of disc-dominated torque is expected as long 
as both M(Ri n ) and Q are roughly constant. 

After the first constant phase, there are three types of 
discontinuity that can occur on the braking index. The most 
dramatic is the transition between spin-down and spin-up 
phases, which formally results in an infinite value for n. 
These reversals correspond to switches between NS surface 
accretion and propeller phases. In our simulations, they hap- 
pen only at the high end of the initial fallback rate (Fig. [2l 
left-hand panels); for lower inflow rates, the system remains 
in the propeller phase all the time. Another discontinuity is 
caused by the transition of the accretion rate at the inner 
disc radius from the early-time Eddington rate to a power- 
law decay. This discontinuity in Cl will appear only if the 
disc torque dominates over the dipole torque at the time 
when M undergoes the transition. Last, as discussed earlier, 
if, at any point during the pulsar evolution, i? m becomes 
formally larger than i?i c , the magnetic pressure of the pul- 
sar is not able to push the disc further out, and hence the 
inner disc boundary tracks R\ c . Since the time dependence 
of Rm (modulated by M(Ri n )) and that of R\ c (dictated by 
Q) are different, these two radii increase at different rates 
with time; hence the inner radius of the disc, Ri n , can tran- 
sition between R m and Ri c (or viceversa) during the coupled 
NS/disc evolution. Any time there is such a transition, the 
torque in the right side Eq.([TJ changes functional form with 
time, and hence so does Q. 

We need to remark that our semi-analytical treatment 
of the torque (Eq. [TJ) naturally leads to the transitions be- 
tween the various phases of the coupled disc/NS evolution 
to be abrupt. In reality, it is likely that at least some of 
these transitions will be smoother. However, we also note 
that sharp transitions in the period (and hence discontinu- 
ities in its derivative) are often observed in bynary systems, 
where the NS dynamics is known to be driven by the disc 
(see Bildsten et al. 1997 for a comprehensive study). 

An important feature in the evolution of the braking 
index is a period of time during which n ~ 1. This phase 
sets in when the disc torque begins to dominate over the 
pulsar torque, and at the same time Q Qk- Then, the 
NS spin is driven by the equation IQ ~ —2MR^ n Q, and 
the braking index can be written as n ~ 1 + (t/t t ) l ' 3 ^ 7a ~ 1 \ 
where t* is a constant dependent on a combination of the 
system parameters. For our choice of a = 19/16, n ~ 1 if 
i > f,, and the two conditions above are satisfied. In our 
simulations, we find that the (short-lived) n ~ 1 phase is a 
tell-tale feature of fast-born pulsars, but is not seen in slow- 
born pulsars, because the latter do not go through a phase 
with Q 3> Qk- 

At late times, we find that the disc torque dominates 
over the dipole torque both for the slow and the fast pulsars; 
Q and Qk track each other, and the torque evolution is ap- 
proximately given by IQ w 2M R^Qk (Rm)[i — Q/Qk (Rm)]- 
The precise time evolution of the torque depends on the in- 
dex a of the accretion rate. For a = 19/16, the time de- 
pendence is close to oc t~ 1,5 , like in the case of a dipole- 
dominated torque, and n m 3. Larger values of a would 
yield somewhat smaller values of n at late times. 
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Figure 3. Time evolution of the various components of the 
NS+disc luminosity, for fast-born (Pq = 5 ms) and slow-born 
(Pq = 300 ms) pulsars, accreting at t = with either a high 
(Mo.max = 10 28 g s" 1 ), or a low (M ,min = 10 25 g s" 1 ) inflow 
rate. L t h indicates the thermal surface emission of the NS. 
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2.3 Luminosity Evolution 

The luminosity of an NS with a fallback disc is given by 
the sum of two terms: a thermal component due to the hot 
surface of the NS, and another thermal component (multi- 
blackbody) due to radiatively-efficient viscous dissipation in 
the accretion flow. The ultimate source of power is inter- 
nal heat for the first term, and the gravitational potential 
drop of the fallback matter, for the second term. In turn, 
the accretion power can be split into separate contributions 
from the accretion disc, from the boundary layer, and from 
the NS surface. Furthermore, when the disc does not pen- 
etrate inside the light cylinder (i.e. when R[ n = Ri c ), the 
pulsar mechanism is likely to operate and then magneto- 
spheric emission, powered by the pulsar spin down, is also 
present. 

More in detail, the various emission components are as 
follows. 

(i) NS emission. The main component, always present with 
or without a disc, is the surface emission of the NS, due 
to its internal cooling. We model the resulting X-ray lumi- 
nosity using the standard cooling model of Blaschke et al. 
(2004), for a 1.4M Q NS, with radius of 10 km. Furthermore, 
as discussed above, there can be periods during the coupled 
NS/disc evolution when the pulsar mechanism is able to op- 
erate, converting a fraction of the spin power of the NS into 
non-thermal X-ray radiation. Observationally, a correlation 
between the spin power E Iot and the 2-10 keV luminosity 
has been investigated in numerous studies (e.g. Saito et al. 
1997; Possenti et al. 2002; Kargaltsev & Pavolov 2008). In 
particular, Possenti et al. (2002) derive the best-fitting rela- 
tion 



logLx,[2-io] = 1-34 log E Iot - 15.34, 



(11) 



with Lx and E Iot = IQ,Cl in CGS units. When present, such 
a component would manifest itself as a power-law spectrum 
of photon index F ~ 2. 

In this work, since we are specifically interested in the 
properties of fallback discs, and their effects on NSs, we com- 
pute in detail the various components of the emission due 
to the disc itself (viscous dissipation, irradiation) and its in- 
teraction with the NS (boundary layer, accretion). However, 
since the surface emission of the NS cannot be hidden, we 
also include this contribution for reference and comparison 
to the components above. For the magnetospheric compo- 
nent, on the other hand, we will limit ourselves to a discus- 
sion, since the presence of this component depends on the 
location of Ri n at a given time, and its magnitude is deter- 
mined by the timing parameters of the NS at that time. 

(ii) Surface accretion luminosity. This term is given by the 
potential and kinetic energy released by matter falling from 
the magnetospheric radius to the surface of the NS; hence, 
in the Newtonian approximation it is: 



Lace = GMM 



Rns 



Rns) 



(12) 



For the characteristic radius of a NS (~ 10 km), the 
effective blackbody temperature of the surface accretion 
luminosity spans the soft X-ray band, between ~ 0.1 keV 
and ~ 2 keV, for a large range of luminosities, from those 
characteristic of Bondi accretion from the instellar medium 



l|Zampieri et alj 1 19951 ) up to the Eddington luminosity. 
It is important to notice that this component is present 
only when R ln = i? m < R co , where R co is the corotation 
radius, i.e. the radius at which the Keplerian rotation of 
the material is equal to the rotation rate of the star Q 
(condition equivalent to < Q^^Rm)). Otherwise, the 
pulsar is found in the propeller phase, accretion is inhibited, 
and the corresponding luminosity contribution can be 
considered negligible. 

(iii) Boundary layer luminosity. This contribution is due to 
the release of energy in the boundary layer between the Kep- 
lerian disc and the magnetosphere. We estimate the bound- 
ary layer luminosity with the following argument, following 
Perna et al. (2006). We adopt the "elasticity parameter" /3, 
which is a measure of how efficiently the kinetic energy of the 
NS is converted into kinetic energy of ejected matter through 
the magnetosphere-disc interaction (Eksi et al. 2005) . In the 
limit of a completely anelastic propeller (/3 = 0), the mag- 
netosphere forces matter to corotate with it during both the 
accretion and the propeller regime. On the other hand, in 
the limit of a completely elastic propeller (/3 — 1), matter 
is slowed down in the boundary layer; this term is present 
only during the propeller phase. 

For a generic value for the elasticity parameter, the lu- 
minosity of the boundary layer can be written as: 



-d 

4tt 



for Rm < R co 
/?) [R 2 m (fl 2 - for R m > R co 



(13) 

Given the uncertainties in the parameter j3, here we adopt 
P = 1, which yields a lower limit to the attainable boundary 
layer luminosity during the lifetime of the NS. 

(iv) Disc luminosity. This contribution is clearly present in 
both the accretion and the propeller reg ime. For an optically 
thick , geometrically th in disc model (Sh akura fc Sunvaevl 
1 19731 : iFrank et all 120021 '). the luminosity in the [^1,^2] fre- 
quency band is given by an integration over radial annuli, 



V2 , hv z 
dv—- 

c A 



RdR 



e hu/kT(R) _ 1 



(14) 



' "1 "J R in 

each producing a blackbody spectrum at a temperature 



T(R) = 



( 3GM NS M(R) \ 
V 8nR 3 a 



1/4 



Bin 



1/2 



1/4 



(15) 



Besides emission due to internal energy dissipation, the disc 
can brighten up also as a result of reradiation of X-rays 
impinging from the central object. We use Vrtilek et al.'s 
(1990) analytical expression for the effective temperature of 
an irradiated disc, under the assumption that the disc height 
h oc r" and irradiation is the dominant form of heating. If 
particles in the disc are in Keplerian motion, then n = 9/7, 
and the temperature profile is given by 



Tx(r) 



2.32 x 10 4 



\Jk B / 'iirriYiLxoJ 



UiraGM 



J_ 
0.5 



2/7 



2/7 



2/7 



r 



3/7 



K, (16) 
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where Lx is the X-ray luminosity of the central source, and 
/ parameterizes the uncertainty in the disc structure and 
albedo. To a first approximation, the total luminosity in the 
[fi, V2} band is the sum of the viscous dissipation and rera- 
diation components; the former dominates at small radii, 
the latter dominates in the outer disc. 

Recall that we are not assuming a constant M through- 
out the disc. As we discussed in Section 2.1, at earlier times 
M(R) ~ R because of wind losses. Therefore, during the 
early super-critical accretion phase, the temperature profile 
in the inner disc is T(R) ~ R' 1 ^ 2 , flatter than the standard 
disc-blackbody profile, which scales as R~ 3 ^ A . Our choice of 
accepting an Eddington-limited accretion rate at the inner- 
most disc radius and a much higher mass inflow rate at the 
outer boundary implies that the integrated disc luminosity 
may exceed Z/Edd at early times. However, it is easy to see 
that for M(R) ~ R, this excess is only a logarithmic func- 
tion of M(R ou t) / Msdd- This is a restatement of the well- 
known result that L disc ~ [1 + ln(M/MEdd)1 during super- 
Eddington accretion (|Shakura fc Sunvae vl 19731 ; iFrank et al.l 
2002). In our case, the highest possible inflow rate at the 
outer radius considered in our model, M(R ou t)/MEdd ~ 
10 10 , would lead to a bolometric disc luminosity ~ 20 Z/Edd, 
at least for the initial transient phase (or less, if R m 3> 10 
km). We accept this possibility of transient super-Eddington 
luminosity as we are modelling a situation characterized by 
massive outflows rather than a steady-state solution. A more 
complex modelling of the accretion flow in that situation 
with a slim-disc rather than a standard-disc plus outflow 
approximation is beyond the scope of this work. 

We begin by exploring the relative magnitude of the 
various components of the luminosity as a function of NS 
age, for a few representative combinations of NS and disc 
parameters. In particular, for the same choice of initial NS 
spin and disc accretion rate as in Fig. [2l we consider the 
time-evolution of the surface accretion, boundary layer, and 
the disc luminosity in the 0.1-10 keV, the 3-4/im, and in the 
80-120 GHz bands. In the IR and the mm band we com- 
pute only the disc luminosity, since the pulsar contribution 
in these bands is not well known. For comparison, we also 
show the timing evolution of the NS thermal emission in 0.1- 
10 keV. We assume an outer disc radiufl i? ou t = 1000 i?NS, 
and a face-on geometry. Hence the predicted luminosities 
represent upper limits for discs with generic inclinations 
with respect to the observer line of sight. 

Fig. [3] shows the evolutionary path of the vari- 
ous components of the luminosity, for the four com- 
binations {Mo.max.Po 7 }' {M),max,fo}, {M ,min , P/ } and 
{Mo,mm,Po}, Where M , m ax = 10 28 g S" 1 and M , m in = 

10 25 g s -1 . We suggest the following physical interpretation 
of these results. At the low end of the initial fallback rate, the 
NS will never be found in an accretion phase; hence, both the 
surface accretion and the boundary layer luminosities will al- 

2 In reality, the outer radius of the disc is expected to be an 
increasing function of time. Modeling this in detail is however 
beyond the scope of our paper, since the magnitude of the disc 
torque depends only on the location of the inner boundary of 
the disc. The X-ray and IR luminosity of the disc is also not 
very sensitive to the outer radius (unless very small), since it is 
produced relatively close to _R m , while the mm luminosity is more 
sensitive to the exact value of -Rout, see e.g. Posselt et al. (2010). 



ways be negligible. In the 0.1-10 keV band, it is interesting to 
note that the disc luminosity exceeds that from the NS sur- 
face under all choices of parameters, for an age t < 10 3 yr. 
The higher the accretion rate, the longer is the time over 
which the disc luminosity exceeds that of the NS (compare 
panels labeled with Mo, max with those labeled Mo.min in 
Fig[3]). On the other hand, the magnetospheric component 
of the luminosity due to rotational energy losses, when the 
pulsar mechanism is active, would produce very high lumi- 
nosities at early times for millisecond pulsars (see e.g. Perna 
& Stella 2004). However, spectral modeling can in principle 
disentangle the thermal disc component from the powerlaw 
emission from the magnetosphere. 

When the initial fallback rate is very high, we have al- 
ready noted (Section 2.2) that, at some epochs, Q,k > fi, 
and the NS is temporarily spun up while accreting. During 
these transient bright phases, the X-ray brightness of the 
system is dominated by the accretion luminosity onto the 
surface of the NS. The exact duration of this transient bright 
phase, and the time at which it sets in, are a function of P, 
M and B; a more detailed discussion is beyond the scope 
of this work. The physical significance of these accretion- 
powered phases is that, if fallback discs around pulsars are 
ubiquitous, and they set in at high accretion rates, then they 
would give rise to a new class of X-ray transients. Further- 
more, we note that a sudden increase in brightness could 
also occur if, at some point during the evolution, the inner 
radius switches from R[ n < Ri c to Ri n = R\ c , at which point 
the pulsar mechanism can turn on again. 

As already discussed in previous work on fallback discs, 
the disc is expected to be brighter at wavelengths longer 
than the X-rays. This is clearly seen by the mid-IR disc lumi- 
nosity (computed in the 3-4 /im band) also diplayed in Fig. [3] 
the mid-IR luminosity is, especially at later times (when 
the disc is cool er), substantially larger than the 2-10 keV 
luminosity fsee lPerna fc HernquistlboOOl ; IPerna et aTll2000l 
for detailed spectra of fallback discs). Searches for fallback 
discs in the optical and IR have been numerous, but they 
have mostly yielded upper limits on the emission. Mignani 
et al. (2007a,b) derived an upper limit on the K luminosity 
of the pulsar PSR J1119-6127 of 6.6 x 10 s0 erg s _1 . A ded- 
icated search for fallback discs in four supernova remnants 
was performed by Wang et al. (2007a). Their IR limits are 
on the order of a few x 10 29 erg s _1 for objects of a few 
kyr of age. Taken at face value, these limits, once compared 
with the results of Fig. [31 would imply that fallback discs, 
if ubiquitous around young NSs, are unlikely to be accret- 
ing at very high rates and/or be very large (recall that here 
we have assumed Rout = 1000 Rns)- Detection of a fallback 
disc around an isolated NS was reported by Wang et al. 
(2006). The mid-IR luminosity was 8.5 x 10 31 erg s" 1 (for 
a distance of 5 kpc). This value is compatible with our pre- 
dictions. More generally, IR emission has been detected in 
a number of isolated NSs (see e.g. Mignani 2011 for a sum- 
mary). However, lack of a detailed spectral characterization 
makes it typically difficult to distinguish between magneto- 
spheric and disc emission. 

Longer wavelength observations (in the submillimeter) 
have been reported by Posselt et al. (2010) for the iso- 
lated NS RX J1856. 5-3754. Using the bolometer array on 
the APEX telescope, they determined a flux limit of 5 mjy 
at 345 GHz. At the estimated distance of 167 pc, this yields 
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Figure 4. Left: the evolution of a representative, short-period NS spin distribution at birth (solid line) for the case in which the NSs 
evolve by dipole radiation alone (bottom panel), and for the case in which they evolve as the result of dipole losses together with 
magnetospheric interaction with a fallback disc (upper panel). Right: same as in the left-hand panels, but including only the pulsars 
above the death line. 



a limit of 5.8 x 10 erg s~ on the luminosity. This value 
is unconstraining for the discs considered here. And in fact, 
Posselt et al. (2010) found that it could only set limits on 
a hypothetical very large outer radius of the disc. The best 
prospect for detecting the presence of discs around young 
NSs remains in the optical/IR wavelength range. 



3 POPULATION PROPERTIES 
3.1 Timing Properties 

We proceed now to simulate how the ubiquitous presence of 
fallback discs affect the statistical timing properties of the 
NS population, for a certain birth spin period and magnetic 
field distribution. The magnetic field B, and the initial disc 
inflow rate Mq are randomly generated according to their re- 
spective distributions discussed in Section 2.2. For the initial 
spin, we consider two distributions: one centred at Po = 5 
ms (representative of fast-born pulsars) and one at Po = 300 
ms (slow-born pulsars); we assume o"i O gP — 0.2 for both dis- 
tributions. Finally, for plotting purposes, we normalize the 
results to a pulsar birth rate N = 1/100 yr -1 . 

In order to simulate the distribution of observable radio 
pulsars, we need to take into account that the radio emission 
fades with age. Here, we adopt the death probability band 



empirically denned by I Arzoumanian et alj (|2002j), centered 
around the line P 15 /P 3 - 10 DL ° = 0, where P is in the sec- 
onds, Pi 5 = P/10 -15 s/s, and DLo = 0.5. The radio emission 
can also be disrupted or modified when the disc penetrates 
the magnetosphere inside the light cylinder (e.g. Shannon 
& Cordes 2009 for an extensive discussion). This does not 
affect the long-term evolution of the period and timing prop- 
erties computed here. However, if disc penetration inside the 
magnetosphere were to dim or completely shut down radio 
emission, then a fraction of pulsars in the sample may not 
be visible. Previous investigations (e.g. Alpar et al. 2001; 
Menou et al 2001a) assumed the inner edge of the disc to be 
always located at the light cylinder, throughout the entire 
evolution of the pulsar, so that radio emission would always 
be ensured. However, when the accretion rate is very high, 
and the ram pressure of the material is larger than the mag- 
netic pressure of the pulsar, there is no compelling reason 
for imposing that the inner edge of the disc does not pene- 
trate inside the light cylinder. Hence, in this investigation, 
we let the disc torque operate at its physically expected lo- 
cation, that is the smallest between P m and Ri c , with the 
understanding that, if P m < Ri c the radio signal might be 
disrupted. With this caveat in mind, in the rest of the paper 
we conform to the standard language in the literature, and 
refer to 'dead' pulsars only as to the fraction of pulsars that 
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Figure 5. Same as in Fig. [4] but for the long-period spin distribution at birth. 



have passed the death line. For a given distribution of ini- 
tial parameters, we calculate period distributions for both 
the full NS population (whether observable or not) and for 
the sub-population of observable radio pulsars (before the 
death- line cutoff). 

First, we examine the evolution of pulsar periods for 
the initial short-period distribution and no death-line cut- 
off. We show (Fig. 01 left-hand panels) period distributions 
that include all the pulsars younger than a certain age, for 
three different ages. The top panel shows the evolution due 
to the combined effect of disc and dipole radiation loss, while 
the bottom panel shows, for the same initial parameters, the 
standard evolution due to dipole losses alone. An interesting 
and perhaps surprising feature of the evolutionary pattern 
of the spin distribution in the dipole+disc case is that it 
becomes bimodal. An inspection of the evolutionary tracks 
of individual pulsars belonging to each peak shows that the 
shorter-period peak is produced mostly by NSs in which the 
dipole torque dominates over the disc torque; hence, the tim- 
ing evolution of those pulsars is only marginally affected by 
the presence of the disc, at least when the NSs are still rela- 
tively young. This can also be seen by comparison with the 
lower panel, in which we plot the expected evolution in the 
absence of disc torques. Another (relatively small) contribu- 
tion to the shorter-period peak comes from NSs in the disc- 
driven spin-up phase, which can bring their spin close to the 
initial value. The longer-period peak of the spin distribution, 
on the other hand, consists of NSs in which the disc torque 



dominates. We have seen (Section 2.2) that disc torques are 
almost always negative: therefore, disc-dominated NSs will 
more quickly reach longer periods. A similar period bimodal- 
ity due to disc torques is also apparent in our simulated dis- 
tribution of observable radio pulsars, before the death line 
(Fig. [31 right-hand panels) . 

The observed period distribution of isolated pulsars 
(ATNF pulsar catalogue: Manchester at al. 2005) is char- 
acterized by a dominant broad peak centred around a pe- 
riod of several hundreds milliseconds, and another, much 
smaller peak, centered around a few milliseconds. If fall- 
back discs are common around isolated NSs, then our re- 
sults suggest that the millisecond peak might correspond to 
the small fraction of the NS population which was born with 
the fastest periods and for which the disc torque was not sig- 
nificant. However, there may also be a contamination from 
old, weakly-magnetized recycled pulsars that were spun-up 
to millisecond periods by disc accretion in a binary system, 
but have since lost their companion star (eg, Lorimer et al. 
2004). 

Our finding that disc accretion leads to a bimodal pe- 
riod distribution is particularly interesting in light of obser- 
vations of NSs in Be/X-ray binaries, which accrete material 
at a relatively high rate (from a companion, rather than from 
fallback). These systems display a bimodal period distribu- 
tion, with one peak centered around ~ 10 s and another 
around ~ 300 s, previously attributed to two types of SN 
explosion (Knigge et al 2011). We propose instead that high 



© 2012 RAS, MNRAS 000, [[]-?? 



Fallback discs in neutron stars 11 




1 1 10 2 1 3 10 + 1 5 1 6 10 7 1 1 1 2 1 3 10 4 10 5 1 6 10 7 

Age (yr) Age (yr) 




/-> : 
r n 



' 



1 3 1 4 1 o 5 

Age (yr) 



Figure 6. Left-hand panels: simulated timing age distribution for the fast-born (top) and slow-born (bottom) populations in the presence 
of fallback discs. NSs above the death line are plotted in cyan, and those beyond the death line in black. The birth rate is a function 
of age bin, in order to have the same number of datapoints in every log bin (namely, a rate of 1 yr -1 for 10 yr < Age < 10 2 yr, 1/10 
yr -1 for 10 2 yr < Age < 10 3 yr, etc.). Right-hand panels: same as in the left panels, but for the braking index distribution. NSs that 
are spinning up are plotted in red. 



rates of disc accretion and the competing effects of dipole 
and disc torques may have contributed to this bimodality. 

We then repeated our Monte Carlo simulation for the 
long-period birth distribution. We find that the period dis- 
tribution remains single-peaked at all times (Fig. [SJl, in 
marked contrast to the short-period birth case. The rea- 
son is that when the initial periods are long, the disc torque 
dominates over the entire life of the pulsars, for the greatest 
majority of pulsars. 

Fig. [6] shows the distribution of timing ages and brak- 
ing indices versus real age for the full pulsar distribution; we 
have also distinguished pulsars above and below the death 
line. A variable birth rate (higher for the younger objects) 
has been assumed in order to sample pulsars of different ages 
more uniformly. The fast-born pulsar distribution shows two 
clusterings, one at n « 3 and another at n » 1: the for- 
mer from the younger and older pulsars, and the latter for 
the intermediate phase discussed in Section 2.2. Note how a 
large fraction of pulsars has n < 3, but there can be some 
very large values. The pulsar distribution with slow initial 
spins, on the other hand, is largely concentrated in the range 
1.5 ^ n ?i 3, plus a subset of spinning-up sources found pre- 
dominantly with n < 0. Timing ages are generally smaller 



than the real ages for the older objects, and larger than the 
real ages for the younger NSs. 

Braking indices have been measured for a number of 
pulsars. Interestingly, indices smaller than 3 are rather com- 
mon, and most often with values between 2 and 3 (e.g. Liv- 
ingstone et al. 2005a, b, 2006), with sometimes values smaller 
than 2 (Lyne et al. 1996). This is the natural range of values 
predicted by fallback discs, as we have shown. 

3.2 Luminosity properties 

We have seen in Section 2.1 that the initial inflow rate from 
a fallback disc can be much larger than the critical value 
at large radii, and the emitted bolometric luminosity can 
exceed the Eddington limit by a logarithmic factor. In our 
model, during this initial phase, the mass inflow at large 
radii declines as (t/to)~ 19 ^ 16 , and most of the inflow is lost 
in a wind when M out > MEdd- We have also shown that the 
decay timescale is to ~ 10 2 -10 4 s for typical disc parame- 
ters. Thus, the supercritical phase of accretion has a very 
short duration. The rapid decline in the disc luminosity is 
due to the combined effect of the decline in the accretion 
rate, and the increase of the magnetospheric radius, which 
defines the inner truncation radius of the disc. Note that the 
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disc luminosity in the X-ray band, being produced in the in- 
ner parts of the disc, is especially sensitive not only to the 
magnitude of the accretion rate, but also to the spin period 
of the NS. The reason lies in the fact that the inner radius 
of the disc, Ri n , being determined by the minimum between 
Rm and R\ c , depends on both M (through R m ), and on Q, 
(through R\ c ). The longer- wavelength emission, on the other 
hand, being produced at larger disc radii, is less sensitive to 
the precise location of the inner radius, and hence it has also 
little sensitivity on the spin period of the NS. 
The characteristic peak disc temperature is 
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With our chosen range of initial inflow rates, we find that 
the peak temperature is ~ 10 6 K ~ 0.1 keV for an ini- 
tial period of time lasting ~ 10 2 -10 4 yr, and then declines 
~ t -0 - 55 afterwards. Thus, fallback discs are always much 
cooler than accretion discs in NS X-ray binaries, even when 
they are both near the Eddington luminosity. This is be- 
cause the magnetic field of an old NS in an X-ray binary is 
generally weaker than in a young pulsai[f|, and the disc can 
extend closer to the innermost stable circular orbit or the 
NS surface, reaching peak temperatures ~ 1-2 keV. This is 
also the reason why disc accretion tends to spin down young 
pulsars and spin up old pulsars. In summary, we expect the 
direct luminosity effect of a fallback disc to be a super-soft 
component in addition to the thermal component from the 
NS surface. A noteworthy result of our analysis is that at 
early times, t < 10 3 — 10 4 yr, the 0.1-10 keV luminosity 
from the disc is expected to exceed the luminosity from the 
NS surface for the range of initial accretion rates considered 
here, M(t = 0) > 10 25 g s _1 . This is especially so for mil- 
lisecond pulsars, since at early times Ri n = Ric, and hence 
the inner radius of the disc is smaller for shorter-period NSs. 
Lack of detection of this disc component in young NSs would 
argue against an ubiquitous formation of fallback discs, or 
that the initial accretion rates are lower than what numer- 
ical simulations suggest. In millisecond pulsars, discs could 
be blown away by the pulsars themselves (Eksi et al. 2005). 

Lastly, to conclude our discussion of the luminosity, note 
that, when the non-thermal magnetospheric component is 
enabled, it would likely dominate in the 2-10 keV band at 
least for the fastest pulsars. In this case, our model suggests 
that the X-ray colours and spectral appearance of a young 
pulsar with a fallback disc truncated at the magnetospheric 
radius would be somewhat similar to those observed in the 



3 The issue of magnetic field decay (Goldreich & Reisenegger 
1992) is still not fully settled as far as its quantitative details. 
However, observations (e.g. see summary in Zhang 2007) suggest 
that both Atoll and Z sources (the two classes of low-mass X- 
ray binaries) have weaker B-ficlds; this observational evidence is 
based on the X-ray spectra, the luminosity at which they switch 
from propeller to accretor, and also on the existence of recycled 
millisecond pulsars, which could not be spun up to millisecond 
periods if they had a high magnetic field. 



largest class of ultraluminous X-ray sources (Feng & Soria 
2011 for a review), dominated by a power-law component of 
photon index ~ 2, with a soft excess at kT ~ 0.1-0.2 keV 
and possibly a radiation-driven disc outflow. This is likely 
to be a coincidence due to the relative scaling of accretor 
mass, accretion rate in Eddington units, and characteristic 
inner-disc radius in the two classes of objects. If the thermal 
component of the luminosity in both classes of objects is 
Eddington limited, young pulsars with fallback discs will be 
one or two orders of magnitude less luminous than ultralu- 
minous X-ray sources with cool discs, and there would not 
be possibility of misidentification if the distance is known. 
Another main difference is that the X-ray luminous phase of 
a fallback disc lasts only < 10 4 yr, which implies that such 
sources must be found very close to their natal supernova 
remnant; instead, ultraluminous X-ray sources are thought 
to be active for > 10 6 yr and they are very rarely associated 
with supernova remnants. Moreover, the characteristic inner 
radius of the standard disc in power-law dominated ultralu- 
minous X-ray sources tends to increase with the accretion 
rate (for example, R[ n ~ M if it corresponds to the spheriza- 
tion radius), in the supercritical regime; as a result, the ther- 
mal c omponent becomes cooler but more luminous ()Sorial 
l2007l : iKaiava fc Poutanen|[2009l ). Instead, the inner radius 
of a disc truncated at the magnetospheric boundary will be 
pushed outwards as the accretion rate declines (assuming 
that the NS magnetic field varies more slowly than the ac- 
cretion rate): R m ~ M~ 2/r , implying T in ~ M 13 / 28 ~ i" ' 55 
and Lboi ~ j^a/ia ^ ^-1.52 ^ i a ^ e ^j meSj ^he ^isc becomes 
too faint and cool to be detectable in the X-ray band. Then 
the best prospects for direct detection are at longer wave- 
lengths, especially in the IR. 



4 SUMMARY AND CONCLUSIONS 

Fallback discs have been suggested to be a common outcome 
of supernova explosions. In this paper we have explored some 
of the consequences of their presence for the timing and 
spectral properties of the associated NS population. To this 
purpose, we have modelled the torque and luminosity of fall- 
back discs around young pulsars, with a range of values for 
the initial spin and dipole magnetic field. Relying on the 
results of numerical simulations of massive star collapse, we 
assumed that the initial inflow rate of fallback material at 
large radii is highly supercritical, but most of the inflow is 
blown away in a radiatively driven wind, so that the accre- 
tion rate at the inner edge of the disc is Eddington limited. 

We considered two initial distributions of spin periods: 
one with a Gaussian distribution centred at 5 ms, and one 
centred at 300 ms, with a birth rate proportional to the SFR 
(assumed constant). We evolved the period distribution of 
the two populations over time, under the combined effect of 
dipole and fallback disc torques. The main outcomes of our 
analysis are summarized below. 

(i) We find that, within the wide range of initial param- 
eters explored, the effect of a fallback disc is mostly a spin- 
down. Matter from the disc is propelled out rather than be- 
ing accreted on the NS surface, except for a few short phases 
when the accretion rate is very high. During these phases, 
which can last up to several thousand years, the NS-disc sys- 
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tem brightens up in X-rays to a luminosity that can reach 
the Eddington value, and the NS is spun up to levels that 
can bring the spin close to its initial value. 

(ii) If fallback discs are common, they substantially affect 
our estimates of the initial spin birth of pulsars (which gener- 
ally assume spin down losses due to dipole radiation alone) . 
In particular, the period distribution of fast-born pulsars is 
predicted to evolve towards a bimodal distribution, unlike 
the distribution in the absence of fallback discs, which re- 
mains single-peaked if it started as such. On the other hand, 
for slow-born pulsars, the main long-term effect is simply a 
shift of the distribution to longer periods. 

(iii) Timing ages as measured for the NS / disc system gen- 
erally overestimate real ages for young pulsars, and underes- 
timate it for old pulsars. Braking indices cluster in the range 
1.5 ^ n £ 3 for slow-born pulsars (apart from a short-lived 
spin- up phase with —0.5 ^ n — 2), and —0.5 n ,$ 5 for 
fast-born pulsars. Younger objects tend to have n £j 3. Large 
values of n, while not common, are possible, and associated 
to phase transitions in the pulsar-disc system. 

(iv) In addition to changing the spin evolution of a pulsar 
(and therefore also its rotation-powered X-ray luminosity), 
a fallback disc can be directly detected as a bright source 
in a broad wavelength range. We found that at early times, 
the characteristic inner-disc temperature kT[ n ~ 0.1 keV 
for typical initial values of the magnetic field, B ~ 10 12 
Gauss, and Eddington-limited accretion, and declines as 
Tin ~ M 13/2S ~ t~ - 55 at later epochs. The disc has a low 
temperature compared with typical accretion discs in (older) 
NS X-ray binaries, because it is truncated at the magneto- 
spheric radius. The 0.1-10 keV luminosity generally exceeds 
that of the NS surface at early times, if the initial accretion 
rates are highly superEddington. Lack of detection of such a 
soft X-ray excess in young objects can help put constraints 
on the ubiquity of fallback discs. 
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